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In the hexagonal columnar phase of chiral polymers a bias towards cholesteric twist 
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•'-j ■ state" with twisted bond order. Polymer trajectories in the plane perpendicular to their 

j^ ■ average direction are described by iterated moire maps of remarkable complexity. 
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A notable feature of biological materials is the profusion of long polymer molecules 
with a definite handedness. DNA, polypeptides (such as poly-7-benzyl-glutamate) and 
polysaccharides (such as xanthan) can all be synthesized with a preferred chirality. Long 
polymers in dense solution often crystallize into a hexagonal columnar phase. When the 
polymers are chiral this close packing into a triangular lattice competes with the tendency 
for the polymers to twist macroscopically [1] as in cholesteric liquid crystals. Similar to 
the twist grain boundary phase of chiral smectics [2] , macroscopic chirality can proliferate 
when screw dislocations enter the crystal. Like flux lines in a type II superconductor, dis- 
locations only appear provided the free energy reduction from the chiral couplings exceeds 
the dislocation core energy. If the chirality is weak, a defect free hexagonal columnar phase 
persists, as in the Meissner phase of superconductors. 

In this letter we explore the effect of chirality on the hexagonal columnar phases [3] 
of long polymers in detail. We neglect for simplicity heterogeneity along the polymer 
backbones and work with a two component displacement field perpendicular to the local 
polymer direction. The usual chiral term relevant for cholesteric liquid crystals produces 
the polymer tilt grain boundary phase, similar to the smectic-A* phase [2]. We find, as 
well, an additional term in the free energy which favors the rotation of the bond order 
along the average polymer direction. This term leads to braided polymers with twisting 
describable by a sequence of moire patterns. 

Using continuum elastic theory we estimate the critical values of the chiral couplings 
above which screw dislocations enter. On a more microscopic level, we propose a set of 
lock-in moire textures which should be especially low in elastic energy. These states are 
entangled and self-similar, and lead to tortuous polymer paths of remarkable complexity, 
reminiscent of chaotic dynamical systems. 

When polymer nematics crystallize the areal polymer density in a plane perpendicular 
to the average direction may be approximated as 

P~Po + X]pa(r)exp{-iGa-f} (1) 



where the {Gq} are the six smaUest reciprocal lattice vectors of a triangular lattice. 
We take the average polymer direction here and throughout to be z, with Ga -5 = 0. 
Each plane wave is modulated by a spatially varying magnitude and phase Pair) = 
|pQ(r)|exp{iGa-'u(r)}, where w is a two-dimensional displacement field. In addition, we 
have a nematic order parameter, n ^ z + 6n and, in the case of a triangular lattice of 
polymers, a bond field 9^ which measures the bond-orientational order in the xy-plane, 
modulo 2n/6. Under a global rotation about the x-axis or y-axis by an angle O^ or 9y 
respectively, h ^' h + 9xh ^ z + 9yX — 9xy, i.e. dn ^' dn + 9yX — 9xy. Similarly under a 
global rotation about the ^-axis by 6*2 , 6*6 — ^ 6*6 + 6*^ . Under such rotations, f ^ f — 9xr, 
leading to a position dependent change in the phase of Poi{r). To insure rotational invari- 
ance, derivatives of pa (r) must be accompanied by the fields 9q and dn [4] and we are led 
to the Landau free energy 
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Crystalline order arises for b sufficiently negative so that (pa) j^ ■ Due to the third order 
term, this transition will, in general, be first order. 

To this free energy we add two more pieces. It is convenient, but not essential, to 
imagine that the hexagonal columnar density waves arise from a phase with local nematic 
order as well as six-fold bond order perpendicular to the director axis. This "N-l-6" phase 
has been studied by Toner [5] . Although the present experimental evidence for N-(-6 order 
in polymer nematics is sketchy, such phases seem highly likely in columnar systems, in 
analogy with the hexatic order expected for vortex lines in high temperature superconduc- 
tors [6] . In the absence of chirality we have the usual elastic energies of the nematic field 
and the bond angle field, 
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where the {Ki} are Frank constants, and K'J^ and K^ are hexatic stiffnesses paraUel and 
perpendicular to z. For long chain polymers Ki ^ ^2,^3 [7]. 

Finally, we add those chiral terms which respect the nematic symmetry. If there is 
no preferred direction along the polymers then the free energy must be invariant under 
h — > —h. Additionally, since changes in 9q are measured with respect to the h axis, under 
nematic inversion, 9q -^ —Oe^ and the vector v = V9q changes sign. Under spatial inversion 
n -^ — n, 9q -^ —9q and hence v -^ v. The chiral free energy density JF* contains two 
distinct terms invariant under nematic inversion but which change sign under parity [8], 

JF* == -K2qon-{'Vxh) - K'\qov-h ^ -7Vx x6n- ^'3^9^ (4) 

with 7 = K2qQ and 7' = i^jj^o- The first term often leads to cholesteric twist with period 
271 /qq along a line perpendicular to the plane in which the nematic director lies. The 
second induces twist into the bond-order parameter 9^ with period 27r/go along z. The 
consequences of this additional chiral coupling is the main subject of this paper. The total 
free energy is F = j d^ [JF^ + JFq + JF*]. Note the close similarity between (2)-(4) and 
the Ginsburg-Landau theory of a superconductor in a magnetic field. The fields 5n and 
9q are "gauge fields" minimally coupled to the complex order parameters {paif)} by the 
constraint of rotational invariance. There are two distinct "magnetic" fields 7 and 7'. 

Provided 6 ^ we can set p^ = |po| expliGa ■'u(r)} and minimize F to find dui = d^Ui 
and 6*6 = ^eijdiUj. The resulting elastic free energy is now 



F = j d\ Lul + Kl + K3(9>0' - 7V, X 5n - id,9t 



(5) 



where Uij = \{diUj + djUi), p = l\G\'^\po\'^{A + B) and A = \\G\'^\po\'^{A - B). In the 
columnar crystal the two chiral terms are the same if dzdi = didz- However, in the presence 
of dislocations derivatives do not commute, V^^x^n = eijdidzUj 7^ dzeijdiUj = dz9Q. Here, 
and throughout, i,j--- refer to indices in the xy-plane. Burgers vectors of dislocations in a 
hexagonal columnar phase lie in the xy-plane, and there are three generic types [3] : a screw 
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dislocation, an edge dislocation with tangent along z, and an edge dislocation lying in the 
xy-plane. The latter defect requires aligned polymer ends which we neglect in this paper. 
The remaining dislocations must lie in a plane spanned by their Burgers vector h and z, 
which amounts to choosing dislocation complexions with axy = ctyx, where the dislocation 
density tensor a^i is the density of dislocations with tangents along the 7-direction with 
Burgers vectors pointing in the i direction [6]. 

Proceeding as in [9,6] we introduce a new field w^i = d^Ui away from any disloca- 
tions. The non-commutivity of the derivatives of u is represented by the dislocation density, 
^fiu-ydiyW^i = —a^i. One can solve for w^i in terms of the dislocation density a^i and find 
the equilibrium displacement field in the presence of crystal defects [10]. In terms of the ne- 
matic and bond order field non-commutivity of derivatives means 2dzOQ — 'Vj^xdn = — Tr[Q;]. 
The energy per unit length of a screw dislocation is finite, while that of an edge dislocation 
lying along the z direction diverges logarithmically with system size [11]. The screw dislo- 
cation energy depends on the dimensionless parameter S = (Ka^^/iU^^)^/'^, where ^^ and 
^z are short distance cutoffs in the perpendicular and parallel directions, respectively. We 
estimate 6^1 for polymer crystals which is the condition that the response to chirality 
be similar to that of Type II as opposed to Type I superconductors [2] . In this limit the 
free energy per unit length of a screw dislocation is /s = {fi^^'^K^' 6^)/(27rV2^^), where b 
is the length of the Burgers vector. 

If the chirality is strong and 7 ^ 7' we expect the polymer analogue of the Renn- 
Lubensky twist-grain-boundary state [2]. Each of these tilt grain boundaries (TGB) is 
composed of a parallel array of screw dislocations lying, say, in the xz-plane, pointing 
along the x-axis and uniformly spaced along z with spacing d. As illustrated in Figure 1 
this dislocation texture causes a discrete rotation cj) = ta,n~^ (b/d) in the average polymer 
direction. The spatial integral of V^x5n is non-zero, while the integral of OzOq vanishes. The 
TGB state appears when the chiral coupling 7 exceeds the critical value 7c = fscrew/b. The 
diffraction signature is similar to that of the smectic-A* phase: if the pitch axis is y so that 



the polymers lie, on average, in the xz-p\a,ne, the structure function would consist of two 
Bragg circles with radius 2n/a lying in a Q^Qz-plane and centered on q = (0, ±27r/v^a, 0) 
where ao is the lattice constant and two Bragg spots at Qy = ±47r/v^ao. These features 
should have widths ~ 27r/(i' along Qy where d' is the spacing between tilt grain boundaries. 
If the rotation angle between successive crystalline regions were a rational fraction of 
27r, the Bragg circles would break up into a sequence of Bragg spots spaced around the 
circle. Although a precise determination of d and d' would require more detailed energetic 



calculations, we estimate, following [2], that d ^ d' ^ yO-o/Qo- 

The second chiral coupling 7' has no analogue in chiral smectics. To find a configura- 
tion of screw dislocations which exploits this form of chirality we search for a dislocation 
texture which produces a displacement depending only on z. The only texture which does 
not produce divergent elastic energy is one in which the only non-zero components of a.ji 
are a^n^ = ayy [6,10]. A honeycomb array of screw dislocations, on average, produces this 
dislocation texture. We now find that for a single such grain boundary V±_ x dn vanishes 
far from the boundary while OzOq does not. This sort of grain boundary thus causes a 
net twist of the hexatic order parameter 9q while imposing no net V^ x dn. When 7' ^ 7 
and the chirality is large, screw dislocations penetrate for 7' > 7c = 2/s/6. By choosing 
the rotation angle across the honeycomb dislocation network to produce a high density of 
coincidence lattice sites [12], we produce especially low strain energies across the boundary. 
The superposition of triangular polymer lattices below and above the boundary forms a 
moire pattern. The superposition of many such "moire sheets" along the ^-axis braids the 
polymers with deep minima in the energy at certain lock-in angles. Figure 1 illustrates 
the mapping of polymers across the moire plane, for a rotation angle (p2 ~ 13.2° with a 
particularly dense set of coincidence lattice sites. Polymers in the lower half-space (cir- 
cles) must be connected to the closest available polymer in the upper half-space (crosses) 
to minimize bending energy. Note that the map has a discrete translational symmetry, 
in the sense that any coincidence site could be a center of rotation. Especially simple 



moire maps arise for rotation angles (j)^ = 2 tan~^[v^/3(2n + 1)], n = 1, 2, . . .. It can be 
shown that aU such angles are irrational fractions of 27r [10] so that the structure never re- 
peats upon iteration. Around each coincidence point there are n concentric rings of helical 
polymers The lattice of coincidence points is also a triangular lattice, but with a spacing 
«n = aoA/l + 3(2n + 1)^/2, where ao is the original lattice spacing. The geometrical origin 
of such energetically preferred lock-in angles has no analogue in chiral smectics. The exact 
choice of lock-in angles and spacing between moire planes must again be settled by detailed 
energetic calculations. 

Upon two iterations of the moire map separating three regions of polymer crystal, the 
first coincidence lattice is rotated with respect to the second coincidence lattice by precisely 
the angle of rotation (/)^. Thus the composite coincidence lattice is the "coincidence lattice 
of coincidence lattices", with lattice constant a^/ao. Moire maps iterated p times lead to 
triangular composite coincidence lattices with spacing an{an/cioY~^, i.e. to ever sparser 
lattices of fixed points with intricate fractal structure in between them. Figure 2 shows 
the projected polymer paths for p = 3 and n = 1. In contrast to the TGB state polymers 
for p ^ 1 are highly entangled and wander far from straight line trajectories. The polymer 
configurations and the dislocations leading to them are shown in Figure 3 for n = 1. Near 
the moire planes both 9q and Vx x dn are nonzero. The center polymer is a fixed point 
of all the maps. Any such fixed polymer has a halo of others twisting around it. In this 
special tube, the nematic order parameter takes on the texture of a double twist cylinder 
as found in the low-chirality limit of blue phases [13]. The moire state takes advantage 
of both double twist energies and the new chiral coupling 7'. Closely related states are 
possible in the chiral N+6 phases. Textures of vortex lines in Type II superconductors, 
similar to Figure 3, may arise when the Abrikosov fiux lattice is subjected to a strong 
supercurrent parallel to the field direction [10]. 

In the inset to Figure 2 we show a random selection of polymer paths projected onto 
the xy-plane representing 99 iterations of the moire map. Though these polymers are 
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clearly influenced by the exceptional center point which is a fixed point of all 99 mappings, 
there are still paths of great complexity such that polymers wander far from their original 
positions. In a defect-free hexagonal columnar phase the most intense Bragg spots appear 
in the Qz = plane at the vertices of a hexagon with radius G = 47r/v3flo- The braided 
crystal described by the iterated moire map should exhibit a Bragg ring at \qj_\ = G in 
the Qz = plane. The ring should be very sharp along q±, and of width ~ 2Tv/d' along Qz, 



where d' is the spacing between moire planes. We estimate that d' ^ an ^ v ao/o'o- As in 
the Renn-Lubensky state [2], rational moire rotations with angle (p = 27rs/t, where s and t 
are relatively prime integers, would cause the Bragg rings to break up into t spots spaced 
out around the ring, with additional structure along Qz- 

The phase diagram in Figure 1 summarizes our conclusions. The hexagonal columnar 
phase is like the Meissner phase of Type II superconductors. The two chiral couplings 7 and 
7' cause screw dislocations to penetrate the crystal above critical strengths 7c and 7^, in 
much the same way as vortices enter Type II superconductors above the lower critical field 
Hci- The TGB phase predicted here for chiral polymers is similar to that already observed 
experimentally in chiral smectics [2]. The braided moire state is qualitatively new, and 
its experimental observation would be of considerable interest. Surprising little is known 
about the intricate trajectories produced by iterated moire maps shown in Figure 2. These 
could be studied experimentally via neutron diffraction in hexagonal columnar crystals 
with a dilute concentration of deuterated polymer strands. Numerical studies of the local 
fractal dimension and Lyapunov exponent are currently in progress. 
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Figure Captions 

Fig. 1. Phase diagram of a ctiiral polymer crystal. Insets are representative tilt (TGB) 
and moire grain boundaries. Shaded lines are screw dislocations. 

Fig. 2. A projected top view of the moire map with rotation angle (f)i ~ 21.8° iterated 
three times. The inset shows 20 random polymers paths resulting from the same 
map iterated 99 times. 

Fig. 3. The moire state. The thick tubes running in the z direction are polymers, while 
the dark lines are stacked honeycomb arrays of screw dislocations. 
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Figure 2 




Figure 3 



